function [v2,v3]=v2v3(v,rho,ay,my)

sigma1sq = @(s) (2+sin(s)).^2;

sigma1_square_bar = sigmasquarebar(my,v);

phiprime_sigma1_Phi = @(y) 1/v^2*quadgk(@(s) (sigma1sq(s)-sigma1_square_bar).*normpdf(s,my,v),-Inf,y,...
    'AbsTol',1e-12,'RelTol',1e-12);

space = linspace(-10,10,10^4);

%f=@(s) (2+sin(s)).^2.*normpdf(s,m,v);

%disp((sum(f(space))-f(-10)/2-f(10/2))*20/10^6);
avarage = 0;

for ii=1:10^4
    avarage = avarage+phiprime_sigma1_Phi(space(ii));
    disp(ii);
end

avarage = avarage - phiprime_sigma1_Phi(-10)/2-phiprime_sigma1_Phi(10)/2;

avarage = avarage * 20/10^6;

v2 = v/sqrt(2*ay)*2*rho*avarage;
v3 = v2/2;

